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Abstract 

The perturbative expansion of static force and potential is reanalyzed con- 
cerning its practical applicability. A well behaved perturbative prediction is 
given by the integration of the renormalization group equation for the coupling 
«qq(/ i — V r ) = (Cf) r 2 F(r). Since the Lambda-parameter of the Nf = theory 
is known from the perturbative prediction contains no free parameter. It is 
confronted with recent non-perturbative results. For a^0.3 where the truncation 
error of the perturbative expression is naively estimated to be moderate, it is 
really quite accurate and large "non-perturbative terms" are excluded. 
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1. The perturbative expansion in the coupling as is the most important the- 
oretical tool for analyzing strong interaction effects in high energy scattering 
experiments. Since the coupling decreases with increasing momentum transfer, 
perturbation theory becomes accurate in the high energy regime. In practical 
applications, in particular in the determination of the running coupling itself, it 
is tempting to apply the perturbative series already where the coupling is not so 
small. Popular examples are the determinations of as from r-lepton decays and 
the perturbative evolution of deep inelastic scattering structure functions starting 
at a renormalization point below 1 GeV. In order to establish the applicability 
of perturbation theory, it would be very desirable to study such processes sys- 
tematically as a function of the energy, but either they involve a fixed energy (r 
decays, hadronic Z-decays) or the precision is not sufficient over a larger energy 
range (deep inelastic scattering, e + e~ total cross section, Adler function). Our 
main phenomenological test of perturbation theory therefore is the overall consis- 
tency of the determinations of as from different processes (see [0] and references 
therein) . 

Complementary information may be obtained from suitable observables com- 
puted as a function of Euclidean external momenta (or distance) using lattice 
QCD. Here an important limitation is implied by the necessarily finite lattice 
spacing a and the corresponding momentum cutoff 0(1/ a). For QCD with dy- 
namical quarks (and in large volume), one currently reaches a -1 ~ 2GeV. This 
limitation may be overcome by considering a finite size effect as the physical ob- 
servable || which defines a renormalized coupling^. In the Schrodinger functional 
framework, the method has been completely developed and results for the running 
coupling are available f4j|)|,|I] J| , including in particular the Lambda parameter ex- 
pressed in terms of the low energy scale r ~ 0.5 fm |7]] in the theory without 
quarks |IJ (pure Yang- Mills theory). These results refer to the continuum limit, 
reached by a controlled extrapolation from finite a. In this theory the scale depen- 
dence of the Schrodinger functional coupling, «sf(/-0, is in remarkable qualitative 
agreement with perturbation theory for a < 0.3 and for a < 0.2 the 3-loop ex- 
pression describes «sf(/-0 within better than 2%. It is an interesting question, 
whether the coupling in this scheme is a special case or whether this is a more 
"general property of the theory" . 

If one restricts oneself to iVf = 0, also observables in large volume but still 
relatively small distances may be computed employing very large lattices (64 4 ). 
We have done so for the potential between static quarks in the fundamental rep- 
resentation, reaching distances of r ~ 0.05fm with small discretization errors and 

1 In finite volume, with no other scale involved but the size of the space-time itself, one may 
of course keep the lattice spacing small compared to this scale. 
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in fact extrapolating to the continuum for r>0.1fm|8|]. In perturbation theory, 
the potential has been computed to two loops P, |iT|p] , |il2l , |i"3f , but at the same 
time the usefulness of perturbation theory even at distances as short as O.lfm has 
been doubted ]ll],[14j . We shall explain below that this is a question of a suitable 



renormalization scheme. When the most natural scheme (defined in terms of the 
force) is adopted, perturbation theory is well behaved at such distances and it is 
interesting to compare perturbation theory to the non-perturbative force. This is 
a rather stringent test of perturbation theory since besides the scale dependence 
of the coupling, its absolute value is predicted by perturbation theory (we remind 
the reader that the Lambda parameter is known). 

A previous exploratory investigation, concentrated on the question, whether 



there are "large non-perturbative" terms in the potential at short distances [15 



as they had been argued to exist [[Tj|. It is not easy to give a definition of "large 
non-perturbative term" . We assume here that roughly the following is meant by 
this statement. 

(i) A certain quantity, here the potential V(r), is considered in a region where 
its perturbative expansion looks well behaved, i.e. the n-loop contribution is 
a small correction and significantly smaller than the (n— l)-loop contribution 
(unless the latter is accidentally small itself). 

(ii) The difference between the full non-perturbative observable and the trun- 
cated perturbative series is much larger than the last term in the series. 

With such a definition, necessarily somewhat phenomenological in character, we 
shall demonstrate below that there are definitely no large non-perturbative terms 
in the potential. To the contrary, perturbation theory works remarkably well 
where the criterion (i) is satisfied. 

In the following, we first investigate the perturbative expressions and find that 
a stable perturbative prediction satisfying (i) in an accessible region of r is given by 
the integration of the renormalization group equation for the coupling a qq (fi). We 
then compare perturbation theory to our numerical results, considering also the 
direct relation between the Schrodinger functional coupling asF and a m (r). For 
completeness we also show the potential itself compared to perturbation theory. 



2. In single scale problems, such as the static potential depending only on the 
separation r, the best perturbative prediction is expected to be the integration 
of the renormalization group equation. This is equivalent to defining a physical 
renormalized coupling, often denoted effective charge fll7l,|l8| . In particular the 
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coupling a qq (/i) may be denned through 



F(r) = — = C F ^p,/*=l/r. (1) 

(It will become clear below, why we here consider the force rather than the po- 
tential.) The running of the coupling defines the associated /3-function, 

^9 = Kg) , 9 = (4™)^ (2) 
with a perturbative expansion 

m S ~ -f{bo + b 1 g 2 + b 2 f + ...} (3) 
bo = j^ (H - |Mr) ,h = ^ (102 - f N t ) . (4) 

The solution of eq. (0), 

A, = tibo?)^™*- 1 '^*** { ~ f° dx [fa + & - fc] }, (5) 

relates the coupling g = g(fi) to the Lambda-parameter. This general expression 
turns into a perturbative one by inserting the expansion eq. ([3]) for the /3-function. 
Truncating in eq. (§) after the term b n -\ and solving eq. ([5]) (numerically) for g 
at given \i (in units of A) defines the n-loop RG solution for the coupling. In 
contrast to the frequently used expansion of F(r) (or V{r)) in terms of a^(fi), 
one does not need to choose the scale \i. For the Schrodinger functional coupling, 
this perturbative prediction (n = 3) has been shown to be rather accurate for 
a < 0.3 by comparison to non-perturbative results [§,[[]] . 

In order to obtain a qq from eq. fl5|) we need to insert the Lambda-parameter 
in this scheme. We start from 

A^|r = 0.602(48) (6) 

referring to the case N{ = and extracted at sufficiently high scale \x where the 
perturbative error is negligible 0. With the coefficient cq (known from PJlOll) in 
the expansion 

ttqqO) = "MS^) + CO^MS^) 2 + Cl^MS^) 3 + • • • ( 7 ) 



we then relate A qq to A-^g via 



Aqq = A^e-/^ . (8) 
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For convenience the ratio of Lambda-parameters is listed in Table [I] together with 
that ratio for other schemes. 

The coupling a qq from 2- and 3-loop RG is illustrated in Fig. [I], using the cen- 
tral value AjyfgTo = 0.602 (the 8% overall uncertainty of this number corresponds 
to a common small horizontal shift of all curves in the figure). The perturbative 
expansion appears quite well behaved up to distances r ~ 0.25fm. At r ~ 0.2fm 
one would expect the 3-loop curve to have an accuracy of about 10%. Since the 
force is completely equivalent to a qq it is given with the same relative accuracy; 
the potential may be obtained by integration of the force. 




Figure 1: Running couplings obtained by integration of the RG with truncation of the 
13- functions at 2- and 3-loop and with = 238MeV @. 



Alternatively one may define couplings ay and ay through the potential 



V(r) = -Cf^^- 2 , u = l/r 
r 



(9) 



and its Fourier transform 



V(Q) = -^C F ^1 



(10) 



We note, however, that the 3-loop coefficients 62 are larger in these cases (see 
Table |1|), in particular in the V scheme. As a consequence, the difference between 



the 2-loop and the 3-loop running coupling in this scheme is only small at very 
short distances and this perturbative expansion appears to be applicable only up 
to a ~ 0.15. This is also illustrated in Fig. [I]. 
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exp(7 — 35/66) 


exp(31/66) 


exp(31/66 + 7) 


0.48811(1) 


6f x (4vr) 3 




1.6524 


2.1287 


4.3353 


0.483(9) 



Table 1: Ratio of A-parameters and 3-loop coefficient of the (3- function for various schemes 
for N { = 0. These results follow from . 

Note that the couplings in the two schemes which are defined in terms of the 
potential are restricted to perturbation theory. Non-perturbatively it is not clear 
how to subtract the self energy term in the potential and, in addition, performing 
the Fourier transformation of numerical data known in a finite range of r is possible 
only in a model dependent way. Also for this reason, a qq is the natural observable 
for the comparison between perturbation theory and non-perturbative QCD. 
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4.1303 


3.0650 


2.5208 




SF 




-0.271(10) 


-1.336(10) 


-1.880(10) 


-4.401(10) 



Table 2: Coefficients / 2 s (s ) for s = A s >/A s and N f = 0. 

3. When two different couplings are known non-perturbatively, it is further of 
interest to study how well perturbation theory predicts their direct relation. This 
means matching the two couplings at finite \x instead of through the Lambda- 
parameter, which corresponds to matching for \i — > 00. The perturbative relation 

a s >(sfJL) = a g (ji) + ff s {s)a s {fi) 2 + ff s {s)a s {^) 3 + ... (11) 
contains a freedom of relative scale, s. Indeed, the choice of s is in general very 



important for the quality of the perturbative prediction 0,513]. The only viable 
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criterion for fixing s appears to be to demand that the coefficients fi(s) are small 
("fastest apparent convergence"). The choice s = so = A5//A5 yields /i(so) = 0, 
and in addition the values of 1/2(^0) | are close to the minimum of 1/2(5) (• The 
coefficients /2(so) connecting selected schemes are listed in Table ||. One observes 
that the SF-scheme is very close to the MS-scheme, the qq scheme is not very far, 
but the other schemes have quite large values of f2{so) in their relation to the 
MS-scheme. In particular, the large coefficient between the MS scheme and the 
V scheme means that the direct expansion of the coordinate space potential in 
terms of a^g (or «sf) is badly behaved, as it has been pointed out in |IT| , |r4| . |i3f . 
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Figure 2: Running coupling in the qq scheme. Results for the continuum limit as well as 
additional points at finite /3, corresponding to finite lattice spacing are shown. In the latter case 
the discretization errors were estimated to be smaller than the size of the symbols. The pertur- 
bative curves use A^gro from |l[ with the dotted lines corresponding to the 1-a uncertainties of 
this combination. 

We emphasize the following point. Although the three different schemes qq, 
V, V differ only by kinematics (differentiation, the Fourier transformation) it 
makes a big difference for the applicability of perturbation theory which one is 
chosen to represent the potential. The analysis of the perturbative series them- 
selves suggests that potential and force should be reconstructed from a qq . 

4. In Fig. ^|we compare the non-perturbative results of [§] to perturbation theory. 
The 3-loop RG expression with A^jg at the upper end of the error bar of eq. (||) is in 
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very close agreement with the non-perturbative coupling. In fact the agreement 
extends up to values of a qq where perturbation theory is not to be trusted a 
priory. For a qq <0.3 our criterion (i) above is satisfied but there is no evidence for 
non-perturbative terms in this region. 
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Figure 3: Test of eq. (|l3|). The uncertainty in the combination /itq has been translated 
into an uncertainty for ^(agp^/so)) and a qq (/z). The non-perturbative values for o:sf(aO are 
constructed from the data of [Q. Errors are smaller than the sizes of the symbols. 



The same conclusion is reached when we consider the relation between agp(/x) 
and a qq (/x) at finite /z: We define the 3-loop expression 



M«sf(/-0) = « SF (//) + 1.336 [a SF (ii)f 



(12) 



such that 



® q M = M"sf(/V s o)) + 0([a S F(/V s o)] 4 ) (13) 

as explained above. Then the non-perturbative values of a qq (/x) and of /i(«sf(/V ^o)) 
are compared in Fig. |3[ If the higher order terms in eq. ( p~3|) are negligible, the two 
different quantities should agree. At a qq (/i) ~ 0.3 a difference is visible but this 
is only about 3 x a 4 , not far from the expected size of the next order term in the 
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series.Q At a qq (/i) ~ 0.2, the difference a qq (/z) — /i(o!sf(/V so)) is not significant 
at all. We conclude that also in eq. ( |T3| ) a large non-perturbative term at short 
distances is excluded. 
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Figure 4: The potential compared to different perturbative expressions. Here, 
r c = 0.54r §. 



Finally we show in Fig. ^ the static potential itself || compared to different 
perturbative approximations. Full line and short dashes are given by 

V(r) = V(0.1hr ) + r dyF(y) , F(r) = C F r 2 a m (l/r) (14) 

JO.lSro 

with the 3-loop and 2-loop RG-solution for a qq . Since we know that the 3-loop 
RG-solution for a qq is accurate, this also hold for V(r) computed through eq. (|H|). 
Again, the full line moves very close to the data points (r < 0.5r ), when Ar at 
the upper end of the error bar of eq. @ is inserted. The long dashes represent 
eq. @ with the 3-loop RG-solution for ay. As it was to be expected due to the 
missing stability of this perturbative expression, it fails in describing the potential. 

2 Note that the next order correction is formally enhanced by a logarithm of a, which 
originates from a resummation of IR divergent terms. It reads (Alog(a) + B)a 4 |^ , |24| , |l"4j . 
While A = 9/(47r) has recently been calculated |25| , p6| , B is not known. The Aa 4 log(a) term 
by itself constitutes a small correction in the figure, which would slightly enlarge the difference 
between /i(«sf(m/ s o)) and a qq (/i). 



8 



A similarly bad perturbative expression (not shown here) is the direct expansion 
of the potential in terms of a^g. 



5. In summary, care has to be taken which perturbative prediction (scheme) is 
adopted to describe the potential. However, perturbation theory does its best in 
the following sense. As usual in an asymptotic expansion, one should first inves- 
tigate the apparent "convergence" by comparing subsequent orders and checking 
that they decrease significantly. If this is not the case, one is obviously outside 
the domain of applicability of perturbation theory or has chosen a bad truncation 
(scheme). According to this criterion the /3-function in the qq-scheme may be 
trusted up to a qq ~ 0.3. Other truncations of perturbation theory for the poten- 
tial that we investigated are applicable for much smaller values of the coupling 
only. Therefore perturbation theory suggests that the qq scheme should be used 
in order to obtain a reliable perturbative expression.^ 

Our comparison with non-perturbative results, obtained in the continuum 
limit of lattice simulations (Nf = 0), does confirm that such a perturbative analysis 
is a good guideline - at least in the case at hand. Of course one should not expect 
miracles when one goes up to a ~ 0.3. At such values of the coupling we only 
confirm that the 3-loop perturbative prediction is good to within about 10% and 
indeed in Fig. |3] one sees explicitly that the truncated perturbative series has 
errors of this order. Similar results have been found for the Nf = 2 coupling in 
the SF-scheme (see Fig. 5 in ||). 

Which lessons can we learn for QCD with quarks? Compared to Nf = 0, the 
relevant perturbative coefficients, 62 and /^(so), which are listed in the appendix, 
are roughly a factor two smaller in magnitude for Nf = 3. This suggests that with 
quarks the perturbative prediction for the potential computed through a qq is also 
applicable up to a ~ 0.3 and furthermore in full QCD the issue of the appro- 
priate scheme is somewhat less important. A direct lattice QCD check of these 
expectations is unfortunately not possible at present and here we had to boldly 
generalize from the Nf = case. In addition, these remarks apply to the mass- 
less theory (we have not investigated mass effects). Current phenomenological 
research concentrates on the application of a velocity dependent potential beyond 
the static limit for phenomenological applications to top-quark physics PTj], On 
the one hand, in this application the potential is needed for quite short distances, 
where perturbation theory is intrinsically more precise on the other hand, 
with the velocity entering as a new scale, this represents a more difficult multi- 

3 Of course, other similarly well behaved truncations of perturbation theory might be found. 
The important point is that a scheme with a large 3-loop coefficient such as V is of no use in 
the region a > 0.15. 
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-5033/18 
247T 2 - 2|9 - ^C(3) 

- 2S*C(3) + fvr 4 - ^fvr 2 

314^2 + 3|9 _ 7|4 C(3) 

-0.275(5) x (4vr) 3 



e 2 

325/54 
f + f C(3) 
f + f C(3) + f vr 2 
+ f vr 2 + f C(3) 
0.0361(4) x (4vr) 3 



e 3 




-^7T 2 
81 

32 _ 

27 81 

-0.00175(1) x (4vr) 



Table 3: Coefficients of eq. (jig) 



scale problem. Indeed the renormalization group has already been applied to deal 
with this complication p9 |. Nevertheless, the lessons learnt in our investigation 
may be useful in this context as well; the type of renormalization group improve- 
ment which we found to increase the reliability of perturbation theory (see Fig. ^) 
has not been applied in [E71 so far. 
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Appendix. In order to ease the comparison of different schemes, we here list 
some perturbative coefficients for general Nf whose numerical values for iVf = 
were quoted in the tables above. Some of them could be taken directly from 
the literature P,p!0|,p!9|,p0|,p4],p!2|,p!3|,p0| , others such as bf 1 had to be computed by 



straight forward algebra. 

The three-loop coefficient of the /3-function can be expressed as 



fr 2 — &2 liVf 



=o + (4^)6 (eiiVf + e 2 Nf + e 3 iV f 3 



with e, listed in Table |3|. 

The one-loop coefficient in eq. (|ll|) assumes the general form 

cS'Sf n \ 



f? 



(ai + a 2 N { )/4ir - 8vr6 log(s)), 



(15) 



(16) 



and vanishes for s = sq = exp ((a± + <22iVf)/(327r 2 &o))- «i and a 2 are listed in 
Table with MS as reference scheme. The other coefficients can be evaluated by 

f 1 ss '(s) = f 1 ss "(s)+fr s \s)^d 



S'St 



bo 



bi 



4-7T 



47T 



(17) 
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S' a\ a 2 

V 31/3 -10/9 

v 31/3 + 227 -10/9-47/3 

qq -35/3 + 22 7 2/9 - 4 7 /3 

SF -1.25562 x (4?r) -0.03986 x (4tt) 

Table 4: Coefficients a, of eq. @ for S = MS 



which reduces to 



fS"5/ e N _ (4tt) 2 r. c/ 



/f S (*o) 



60 



6f - 6? 



when we set s = s . 
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